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Conformal Surface Morphing with Applications 

on Facial Expressions 

Mei-Heng Yueh, Xianfeng David Gu, Wen-Wei Lin, Chin-Tien Wu, and Shing-Tung Yau 


Abstract —Morphing is the process of changing one figure into another. Some numericai methods of 3D surface morphing 
by deformabie modeiing and conformai mapping are shown in this study, it is weii known that there exists a unique Riemann 
conformai mapping from a simpiy connected surface into a unit disk by the Riemann mapping theorem. The diiation and reiative 
orientations of the 3D surfaces can be iinked through the Mobius transformation due to the conformai characteristic of the 
Riemann mapping. On the other hand, a 3D surface deformabie modei can be buiit via various approaches such as mutuai 
parameterization from direct interpoiation or surface matching using iandmarks. in this paper, we take the advantage of the unique 
representation of 3D surfaces by the mean curvatures and the conformai factors associated with the Riemann mapping. By 
registering the iandmarks on the conformai parametric domains, the correspondence of the mean curvatures and the conformai 
factors for each surfaces can be obtained. As a resuit, we can construct the 3D deformation fieid from the surface reconstruction 
aigorithm proposed by Gu and Yau. Furthermore, by composition of the M5bius transformation and the 3D deformation fieid, the 
morphing sequence can be generated from the mean curvatures and the conformai factors on a unified mesh structure by using 
the cubic spiine homotopy. Severai numericai experiments of the face morphing are presented to demonstrate the robustness of 
our approach. 

Index Terms —conformai mapping, surface morphing, surface matching 
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1 Introduction 


T he metamorphosis between two objects is com¬ 
monly called "morphing". It is the process of 
changing one figure into another. In recent years, 
image morphing techniques have been widely used 
in the entertainment industry. Many techniques have 
been developed to achieve a desired morphing effect. 
[17] 

In 2D image morphing, S.-Y. Lee et al. [12] proposed 
a technique that generates a -continuous and one- 
to-one deformation from the two-dimensional po¬ 
sitional constraints by using a deformable surface 
model. The transition behavior can be controlled by 
assigning the transition curves for selected points on 
an image. Similarly, A. Gregory et al. [6] proposed 
a method of geometry morphing, which creates a 
morph by defining morphing trajectories between the 
feature pairs and by interpolating them across the 
merged polyhedron. Attributed to these technologies, 
the user can achieve a satisfactory visual effect by 
simply setting the correspondence between graphical 
features. In particular, for images of human faces, 
it is possible to pick the feature points of the face 
automatically. V. Zanella et al. [19] use the so-called 
"Active Shape Models" [3] to find the facial features 
in the 2D images and perform the morphing of face 
images in frontal view automatically. To make the 
metamorphosis look realer, the morphing path must 
be constrained by some physical laws. Researchers, 
such as Verbeek et al. [16], Y. Bao et al. [1] and L. 
Younes et al. [2], proposed methods of determining 
the morphing paths by optimizing an energy func¬ 


tional, which characterizes the intrinsic deformation 
of the surface away from its rest shape. Althrough 
2D image morphing technique has pretty mature, 
3D image morphing remains challenges, especially 
when the virtual real morphing effects are desired. 
In addition, in order to achieve a satisfactory visual 
effect, the texture images also need to be computed 
in the process of visualization. 

On the other hand, with the advance of the three- 
dimensional imaging technology, surface morphing in 
3D has become very important. Comparing to the 2D 
image matching problem, surface matching problem 
is much more difficult, since the surface matching 
involves the correspondence in coordinates and the 
geometric information of images in is far richer 
than images in 2D. D. DeCarlo et al. [4] proposed 
a method of generating smooth-looking transforma¬ 
tions between pairs of surfaces that may differ in 
topology by specifying a sparse control mesh on each 
surface and by associating each face in one control 
mesh with a corresponding face in the other. Y.-S. Liu 
et al. [14] proposed a morphing method that works 
well even when two input 3D triangle meshes have 
very different shapes by creating consistent meshes 
for the original source and target models. Schroder et 
al. [13] proposed a variational approach based on min¬ 
imizing bending and stretching in which correspond¬ 
ing feature points and line segments are matched. E. 
Jeong et al. [11] proposed a feature-based morphing 
technique for two objects equipped with surface light 
fields by using spherical embedding of meshes. M. 
Yang et al. [18] proposes a realistic 3D morphing 
method based on GPU for real-time animation of 
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facial expressions by rendering the real texture on the 
digital character. 

In this study, we propose an efficient way to obtain 
the desired morphing effects. Our goal is to gener¬ 
ate the morphing sequence between human faces in 
which features of the given surfaces are maintained 
during the morphing process. Moreover, in order to 
control digital models and perform digital acting via 
real agents, we also construct the single mesh of 
human faces through the correspondence matching. 
In our approach, we utilize the conformal mapping 
technique. A Mobius transformation and a deforma¬ 
tion from a modified thin-plate model are computed 
to achieve high accurate feature correspondence. A 
single mesh based on geodesic paths among feature 
points is computed very efficiently. With the help 
provided by this single mesh, the texture information 
and the geometric information, such as the conformal 
factor and the mean curvature Hi,i = 
of each frame in the morphing sequence can be 
easily interpolated from the given surfaces through 
various homotopy methods. Then the 3D images of 
each intermediate frame can be computed from the 
reconstruction algorithm proposed by Gu and Yau. 
Our method guarantees the piecewise smoothness of 
homotopy path. So, we can avoid the vibration of each 
frames caused by the unstableness of the 3D camera. 
In addition, we merely require 10% the number of 
original frames to simulate the real motion of human 
facial expressions, so that the size of the information 
data is significantly reduced but the resolution of the 
3D images does not decrease. 

In the following, we review the computation of the 
Riemann conformal mapping in section 2. We present 
our surface matching technique in section 3 and in¬ 
troduce the way of homotopy in section 4. After that, 
we show how we reconstruct the morphing sequence 
in section 6 and demonstrate some morphing results 
in section 7. 

2 Idea for Computing Conformal 
Mapping 

The Riemann conformal mapping plays an important 
role in the surface matching. In this section, we briefly 
review the method of computing Riemann conformal 
mapping. 

2.1 Spherical Conformal Mapping 

The spherical conformal mapping is first computed 
by Gu and Yau [8] in 2003. The idea is based on mini¬ 
mizing the harmonic energy through a nonlinear heat 
diffusion process as following: Suppose the desired 
map is (p : ^ S^. Let ip{v) and n((p(i;)) denote the 

image of the vertex v e A4 and the normal at (p(v), 
respectively. The normal and tangent components of 
A(p are defined as 

(A(^(v))-^ = (A¥3(v),n(</?(v)))n(</?(v)) 



Fig. 1. Spherical conformal mapping 


and 


(A(^(t)))ll = A(^(t)) - (A(p(v))-^, 


respectively. The harmonic map is then computed by 
minimizing the harmonic energy associated with cp 
through the nonlinear heat diffusion process 


d(p 

dt 


-(A(^)ll, 


with the constrain ^ The equation can be 

written in the following form. 


dip 

dt 


-(A^)ll 

-(A(^- (A(^,n((^))n((^)) 
-(A(^- (A(/?,</j) ip), 


and can be solved numerically by using the quasi- 
implicit Euler method (QIEM) [7], 


I + 


{k - 




(m+1) 




(m) 


where, K is the discrete Laplacian, is a diagonal 
matrix with 


K)) K)) • 

The initial map ip^^^ for the above iterative formula 
can be obtained from the Gauss map. It is well known 
that the efficiency of the explicit scheme is usually not 
satisfactory since the time step is generally very small 
due to the diffusive nature. The quasi-implicit Euler 
method [7] proposed by W.-W. Lin et al. is shown 
to be very robust on solving the nonlinear diffusion 
equation. 


2.2 Riemann Conformal Mapping 

The idea for computing Riemann conformal mapping, 
shown in Eigure 2, was first proposed by Gu and Yau 
[10]. Eor a given surface M with single boundary, we 
make it into a closed surface M by using the double 
covering technique. Then we compute the spherical 
conformal mapping of M. Einally, we cut out the 
semi-sphere along the equator and map the unit semi¬ 
sphere onto the unit disk conformally by using the 
stereographic projection. 
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Doubl^overing Stereograph^ Projectiorr 

Mobius Transform 



Fig. 2. Idea for computing Riemann conformal map¬ 
ping 



Fig. 3. The results of the Riemann conformal mapping 


2.3 Numerical Results 

We demonstrate the robustness of the quasi-implicit 
Euler method by computing the Riemann conformal 
mappings of human facial expressions. Two different 
facial expressions and the associated conformal map¬ 
pings are shown in Figure 3. In order to check the 
conformality of the QIEM, we paste the checkerboard 
grid on the image of the Riemann conformal mapping 
and put it back to the surface M. by using the 
inverse of the Riemann conformal mapping If 
the mapping is angle-preserving, every angle should 
be nearly 90 degrees. The histograms of the angle 
distribution, shown in Figure 3, indicate that the 
QIEM is very accurate on angle-preserving. 

3 Surface Matching with Minimal 
Global Error 

Surface matching plays a critical role in surface mor¬ 
phing. L. Younes et al. [2] proposed a simple method 
to interpolate 2D landmark matching by constructing 
a diffeomorphism between two domains. When it 
comes to surface matching, it would be much more 
difficult. However, with the Riemann conformal map¬ 
pings, we reduce the surface matching problem 
into the unit disk matching problem. Hence we can 




Fig. 4. The idea of surface matching 


apply the similar idea of the 2D landmark matching 
to the surface matching. 

For the surface matching, X. Gu et al. [8] use 
landmarks to obtain the optimal Mobius transforma¬ 
tion as the matching function between the spherical 
conformal mappings of two brains. When it comes to 
the matching between two different facial expressions, 
a slight twist is allowed since the transformation of 
facial expressions is actually not conformal. In the 
following, we propose to match the landmarks of each 
facial expressions by composition of the Mobius trans¬ 
formation and deformation from the plate matching. 

We assume that the correspondence of the bound¬ 
ary points of two different human faces Sa and are 
known. The assumption can be realized by putting 
markers on the boundary of the human faces during 
scanning. Suppose m landmarks are selected for each 
facial expression and denoted by 

C Sa 

and 

C Sb. 

Our purpose is to construct a matching function / : 
Sa Sb, where / ^ i = 1,..., m. First, each 

surface is mapped to the unit disk D by the Riemann 
conformal mappings 

• Sd —y ID and pb \ Sb —^ ID. 

We denote ^pa {p^^) by p^^^ = (pdIi,Pd* 2 ) and pb (gb)) 

by 90 , 2 )' respectively, i = 1,..., m. Next, 

a uniform grid points of an n x n checkerboard 
C ranged on A = [—1,1] x [—1,1], j = l,2,...,n^, 
is lay on D. In the thin-plate model, the deforma¬ 
tion field is approximated by the span of the Green 
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Fig. 5. The deformation field of /b = /m o fg where fm 
is the optimal Mobius transformation (OMT) and fg is 
the global matching function (GMF) 



OMT 

OMGMF 

EbUb) 

3.8820 X 10“^ 

6.8836 X lO-'" 

Elocif) 

1.5029 X 10-^ 

4.0888 X 10-^ 

E{f) 

4.8951 X 10“^ 

1.9131 X 10-^ 

TABLE 1 


Comparison between the optimal Mobius 
transformation (OMT) and the optimal Mobius 
transformation with the global matching function 
(OMGMF) 


functions log r of the bending operator at each 
grid point where r is the distance between and 
X e A. Therefore the matching function between the 
unit disks /b : ID) ^ B is defined by /b(^i,^ 2 ) = 
(/b^^ (xi,X 2 ) (^ 1 ,^ 2 )) with X = {xi,X2) and 



2 

log 


X — c 


U) 


for k = 1,2, where are unknown coefficients, j = 
1 ,..., To determine these coefficients for matching 
landmarks on conformal parametric domain, we solve 
the least square problems 


arg min 

a(^) 


— gB,/c , 


k = 1,2, 


where 






log 


9 !'^ - 


i = l,...,m, j = 1, 


Aa ^ and are the conformal factors, resulted from 
the Riemann conformal mappings Lpa and (p^, at 
respectively, and 


a 






QB,k = 


(1) (m) 

B,/e’ • • • ’^B,/c 


The least square problem can be easily solved by 
QR method when < m or solved by Tikhonov 
regularization [5] 

{si + S) qB,k 

when ri? > m. Then the matching function between 
Sa and Sh can be obtained by taking f = o 
Figure 5 shows the result of global matching and 
Table 1 shows the comparison between the opti¬ 
mal Mobius transformation (OMT) and the optimal 
Mobius transformation with the global matching func¬ 
tion (OMGMF). 

There is a variety ways to measure the distortion 
of the matching function, such as Euclidean norm on 
the space, the infinite norm on the unit disk, etc. In 
order to verify that the matching pairs are mapped 
to by /b, i = 1 ,... ,m, we define the conformal 
matching energy 


^b(/b) 




Ai) 

i=l 




(0 (d ^ ■ 

^^B ’ ^B ^ ■ 


which measures the square of Euclidean distance at 
each matching pair on the conformal disk. On another 
point of view, we define the local matching energy 


m 


e 5a, g(') e St, 


which directly measures the square of Euclidean dis¬ 
tance at each matching pair in In order to ensure 
that the global error is relatively small, we also com¬ 
pare the global matching energy, which is defined by 

E{f)= f \f {Sa{x)) - Sb{x)f dx. 

JB 

According to Table 1, the global matching function 
significantly reduces the matching energies, which 
indicates that the global matching function works well 
on the surface matching. 


4 Surface Morphing 

The effect of the traditional image morphing by using 
the direct interpolation is not satisfactory since there 
are shadows in the area of wrong correspondence. In 
3D morphing, it could be even worse. To improve this 
phenomenon, D. Smythe proposed the mesh warping 
technique [17], which partitions the 2D image into 
several pieces and interpolate them piece by piece. 
Base on this idea, we compute a geodesic frame M 
on the surfaces in 3D. The frame M consists of the 
discrete surface geodesics that connect preselected 
feature points. The frame M can be refined to better 
capture characters of various facial expressions. The 
geodesic frames for the eight facial expressions are 
shown in Figure 6. 

The method of computing discrete surface 
geodesics on triangular meshes was proposed 
by D. Martinez et al. [15] which computes the initial 
path by using Sethian's fast marching method and 
correct the path with the path correcting iterations. 
For a mesh with a large number of vertices, the fast 
marching method could be very time consuming. To 
improve the efficiency, we calculate the initial path 
by 

(i) construct the frame on the unit disk, where fea¬ 
ture points are connected by straight line seg¬ 
ments. 
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neutral smile 


bitter smile pain 

Fig. 6. The geodesics on a 
facial expressions 



sad pout 


wry ferocious 

human face with different 







bitter smile pain wry ferocious 

Fig. 7. The partition mesh of the unit disk for different 
facial expressions 


(ii) the initial paths are obtained by taking the in¬ 
verse conformal map of these line segments. 

(iii) Apply Martinez's algorithm to obtain the 
geodesic frame. 

The resulted geodesic frame is called the single mesh. 
In our calculation, the initial paths in the frame 
mostly converge to the geodesics within 5 steps in 
the path correcting iterations. The geodesic frame 
is remapped to partition the conformal parametric 
domain as shown in Figure 7. 

5 Surface Registration via the 

Geodesic Frame of the Surface 

To build an one-to-one surface registration, we rely on 
the aforementioned partition mesh Let (pa{Ma) 

be the partition mesh on (pa{Sa)- The partition mesh 
on (pb{Sb) is determined by the / o (pa(Ma) where / 
is the matching function discussed in section 3. To 
introduce our registration map, first let us introduce 
some notations. The partition mesh of the unit disk is 
denoted by Lp{M){V,F), here 



are the set of vertices and the set of triangles in ip{M) 
where , j = 1,2, is the coordinate of the i-th vertex 
and j = 1,2,3, is the indices of the vertices of 
the i-th triangle, #(V) and denote the number 

of vertices in V and the number of triangles in T. 
Obviously, the closed region determined by Fi can be 
easily represented by 


ai{v)V n) Fa 2 {v)V { 2 ) + n 3 (i;)l/ o), 

^ i ^ i ^ i 

3 

'^ai{v) = 1, ai{v) >0,i = 1,2,3, 

i=l 

where {ai{v), a 2 {v), a^{v)) is the barycentric coordi¬ 
nate of V with respect to the vertices V^( 2 ) 

and V^(3). A simple piecewise linear registration map 
: Sa ^ Sb can be easily constructed by the affine 
mapping 


^(i;) = ai{v)f{V^a))Fa2{v)f{V^,2,)^as{v)f{V^,s,), 

^ i i ^ i 

if G Fi, between each pair of triangles Fi G (Pa(Ala) 
and ~]{Fi) G / o ipa{Ma)r where 


= 2 
a2(v) = i 

asiv) = i 






vV^(3) X vV^(i) 


X vVj,(2) 


Recall (i7. A) is a unique representation of S. So, the 
surface registration can be realized through 




Similarly, the texture images Ta and % of the surface 
Sa and Sb, respectively, can be registrated by 

^^[ra{v)] = rbmv)). 


In the following, we introduce how we utilize 
the above surface registration method to generate 
the morphing sequence through the cubic spline 
homotopy of the mean curvatures and the confor¬ 
mal factors. Suppose 3D images of facial expressions 
So, Si, Sn, are captured a time to, ti,... Using 
the above surface registration method, the registration 
maps : Si-i Si, i = 1,2, ...,At, can be 

easily computed. Using these registration maps, a 
morphing path F{v,t), t G [to,tN] and v G S'o/ can 
be created, here P{v,t) denotes the location where a 
point V e So is morphed at time t. Since {H, A) is a 
unique representation of a surface, the morphing path 
can also be uniquely determined by the evolution of 
the conformal factor and the mean curvature. Here, 
we employee a piecewise cubic spline homotopy to 
interpolate 

{{H, X)o{v), {H, {H, A) 2($2 o 

..., {H, X)m{^n ° ^jv-i o • • • o $i(v))} . 
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A sketch in Figure 8 illustrates this idea. Let 
y[xi,... ,Xn]{t) denote the piecewise cubic spline 
function with given data xi,... The conformal 
factor A and the mean curvature H at F(^, t) can now 
be evaluated by 

{H, X){P{v, t)) = y [{H, A)oM, {H, A)i(^i(^)), 
(F^,A)2(^2 0 ^ 1 ( 1 ;)),..., 

{H, X)n{^n o ^n-1 o • • • o T>i(i;))] (t), 

Similarly, suppose the texture images % of the surface 
Si, i = 1,... ,n, are given. The texture image associ¬ 
ated with the surface along the morphing path P{v, t) 
can be computed by the cubic spline homotopy 

r(P(^,t)) = ^[ro(^),ri(^i(^)),r2(^2 0 ^ 1 ( 1 ;)), 

. . . , Tn{^N o ^N-1 o • • • o (t). 

Finally, by reconstructing the 3D surfaces 
St,t e [to,tAr], from their {H,X) representation 
and applying the computed texture image 7^ to 5't/ 
the morphing sequence between Sq and Sn can 
be obtained for any period from to to In the 
following section, we shall introduce the surface 
reconstrution algorithm in detail. 



Fig. 8. The cubic spline homotopy of the mean curva¬ 
ture and conformal factor of a vertex 


6 Laplace-Beltrami Surface Recon¬ 
struction 


Laplace-Beltrami equations 

AsS{u,v) = 2H{u,v)n{u^v) 

^ ^ ^ =dS 


where 


^ X^{u,v) ^ dv‘^ ) 


and n(it, v) is the normal of the surface S. The detail 
algorithm for solving the Laplace-Beltrami equations 
can be seen in Algorithm 1. 


Algorithm 1 Laplace-Beltrami Surface Reconstruction 
Input: the mean curvature and the conformal factor 
(i7, A) of the surface S and the boundary dS od 
the surface S. 

Output: the surface S. 

1: Set the initial surface normal 
2 : repeat 

3: Solve the boundary value problem 


4; 


where boundary dS is known. 
Update the surface normal 


,(j+i) - 


^(i+i) y c(i+i) 
A2 


5: until convergence 


We reconstruct 7 frames from the video at time 
t = 0.0,0.3,0.6,1.0,1.3,1.6,2.0, respectively. The sur¬ 
face captured from the video at time t is denoted by 
St and the reconstructed surface S{H,X)t is denoted 
by Sf. In order to verify that the reconstructed surface 
is approximately the real surface, we compute the dif¬ 
ference between the real surface and the reconstructed 
surface in sense 




and in sense 


In the previous section, we have obtained the mean 
curvatures and the conformal factors {H,X) at each 
time period U, i = l,2,...,n. Gu and Yau pro¬ 
posed that any surface in 3D Euclidean space can 
be determined by its conformal factor and mean 
curvature uniquely up to rigid motions [9]. We could 
generate the morphing sequence of surfaces {St^ = 
S{H, X)t^}'i^i by reconstructing the unique surface 5't. 
via i = l,2,...,n. 

In the following, we introduce how we utilize 
the mean curvatures and the conformal factors to 
reconstruct surfaces by solving the Laplace-Beltrami 
equations. Applying Gu and Yau's method, for the 
given mean curvature and the conformal factor (i^. A), 
we reconstruct the unique 3D surface S by solving the 


ll-^t - S'tlloo = max 


St{x) - St{x) , 


respectively. Table 2 shows the reconstruction error in 
both 1/2 norm and L^o norm. 


7 Numerical Results of Surface 
Morphing in 3D 

In the following, we show some surface morphing re¬ 
sults via merely two images by using the surface mor¬ 
phing technique which we have mentioned above. In 
each figure, the left most image is the role of initial 
surface Sa while the right most image is the role of 
terminal surface Si,, and the images in the middle are 
the morphing sequence between Sa and Si,. 
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^0 ^ 0.6 ^ 1.3 ^2 



t 

l|5t-5t||2 


0.0 

3.8900 X 10“^ 

6.5598 X 10-* 

0.3 

2.5601 X 10“^ 

3.8562 X 10“^ 

0.6 

4.6325 X 10“^ 

5.1009 X 10-* 

1.0 

2.9394 X 10“^ 

3.9610 X 10-* 

1.3 

4.4131 X 10-^ 

6.5401 X 10-"‘ 

1.6 

2.6350 X 10“^ 

6.6641 X 10-* 

2.0 

4.6418 X 10“^ 

4.4197 X 10-4 


TABLE 2 

Reconstruction error between the real surface St and 
the reconstructed surface St 


Figure 9 and Figure 10 show the morphing se¬ 
quence between two different facial expression, re¬ 
spectively. Figure 11 and Figure 12 show the morph¬ 
ing sequence between faces of two different people, 
respectively. Figure 13 shows the morphing sequence 
between a human face and a loin's head. 

It is interesting that how much the global matching 
function improves the morphing sequence. The affine 
mapping constructed by using the optimal Mobius 
transformation is denoted by and the affine 
mapping constructed by using the optimal Mobius 
transformation with the global matching function is 
denoted by i = 1,2. In order to measure the rate 
of improvement, we reconstruct the surfaces 

sr^s{H,\)r 

and 

where 


t 

1 

to 

l|5f-5d|2 

Improvement Rate 

0.3 

1.4980 X 10-^ 

8.6761 X 10-^ 

42.08% 

0.4 

2.0554 X 10-^ 

1.1229 X 10-^ 

45.37% 

0.5 

2.3328 X 10“^ 

1.0679 X 10“^ 

54.22% 

1.3 

4.1535 X 10-^ 

1.8688 X lO-''^ 

55.01% 

1.4 

3.9139 X 10-^ 

1.9036 X 10-^ 

51.36% 


TABLE 3 

Comparison between the optimal Mobius 
transformation and the optimal Mobius transformation 
with the global matching function 



Fig. 9. The morphing sequence of eye blinking 



Fig. 10. The morphing sequence of mouth opening 


(ff, A)f = ^ [(ff, A)o(^;), {H, A)i($f (r;)), 

o^fiv))] it), 

and 

{H, A)f = ^ [{H, XUv), {H, A)i($f (r;)), 

(if,A)2($?o$f(t,))] it), 

t G [0,2]. Then, we compute the surface difference 



Fig. 11. The morphing sequence from a girl’s face into 
a boy’s face 


\\Sf-Sth = 



and 


l|Sf-S.ll2 = 




- Stix) 



The rate of improvement is defined by 

||gf -5t||2-||gf-5t||2 

\\St^-Sth 


Table 3 indicates that the global matching function 
improves approximately 50% of the surface difference. 

Moreover, extrapolation can also be achieved by 
using this homotopy technique. The real surface at 



Fig. 12. The morphing sequence from a girl’s face into 
another 



Fig. 13. The morphing sequence from Max Planck’s 
face into a lion’s head 































t 

l|St-St||2 

||Sf-St||oo 

1.4 

1.3639 X 10“^ 

2.9274 X 10“^ 

1.5 

1.2816 X 10“^ 

3.0208 X 10“^ 

1.6 

2.6365 X 10-^ 

2.4064 X 10-=^ 

1.7 

3.2057 X 10“^ 

2.3978 X 10“^ 


TABLE 4 

Extrapolation error between the real surface and the 
reconstructed surface 


time t is denoted by St and the reconstructed surface 
at time t is denoted by St = S{H, X)t where 

{H, X)t = y [{H, X)o.e{v), {H, A)i($f (t;)), 
(if,A)i.3($f o$f(^;))] (t). 


Table 4 shows the surface difference between the real 
surface St and the reconstructed surface St = S{H, X)t 
in sense 


||^t-^t|| 2 = f / St{x)-St{x) 
\Jb 


dx 


and in sense 


\\St - S'tiloo = max St{x) - St{x) 


ccGB 


respectively, at time t = 1.4,1.5,1.6 and 1.7. 


8 Conclusion 

In this paper, we proposed a 3D surface morphing 
method between different simply connected surfaces 
with single boundary in which smooth transient on 
both geometric characteristics and texture of the sur¬ 
faces is considered. Similar to the traditional morph¬ 
ing approaches based on boundary representation, a 
wrap has to be created via feature correspondence and 
interpolation between shapes based on the wrap is 
employed to generate the morphing sequence. By tak¬ 
ing advantage of the conformal parameterization and 
the unique surface representation of conformal factor 
and mean curvature, the wrap can be easily obtained 
by the composition of deformations from the Mobius 
transformation and the thin-plate matching function. 
To mimic the non-isomorphic risk that usually occurs 
in matching largely deformed surfaces, a single mesh 
based on geodesic frame is employed. As a result, 
the correspondence, including geometric information 
and texture information, of the whole surface can be 
defined and interpolation among original surface and 
target surface can be computed by the usual cubic 
spline homotopy in a disk parametric domain. Finally, 
the morphing sequence can be generated from the 


surface reconstruction algorithm in section 6. To make 
the proposed morphing approach more attractive in 
real applications, we improve the efficiency in com¬ 
puting the conformal parameterization and geodesic 
frames. We propose an non-linear iterative surface 
reconstruction algorithm (Algorithm 1). The surface 
reconstruction algorithm can be accelerated by using 
the multigrid method on a uniform mesh by which 
multi-resolution surfaces can also be obtained. Several 
morphing effects among different 3D facial expres¬ 
sions are presented to demonstrate the feasibility of 
the proposed morphing method. 
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